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ABSTRACT 

This  paper  describes  the  integrated  modeling  of  a  permanent  magnet  (PM)  motor  used  in  an 
electromechanical  actuator  (EMA).  A  nonlinear,  lumped-element  motor  electric  model  is 
detailed.  The  parameters,  including  nonlinear  inductance,  rotor  flux  linkage,  and  thermal 
resistances,  and  capacitances,  are  tuned  using  FEM  models  of  a  real,  commercial  motor.  The 
field-oriented  control  (FOC)  scheme  and  the  lumped-element  thermal  model  are  also  described. 

INTRODUCTION 

The  development  of  all-electric  aircraft  is  a  high  priority  in  the  avionics  community.  Current 
aircraft  use  a  combination  of  hydraulic,  pneumatic,  and  electric  systems  for  flight  control. 
However,  the  expectation  for  future  airplanes  is  a  single,  electric  system  using  electromechanical 
actuators  (EMAs).  Such  a  system  would  reduce  the  cost  to  build,  operate,  and  maintain  aircraft. 
It  would  also  make  aircraft  lighter,  more  reliable,  safer,  and  more  easily  reconfigurable,  reducing 
the  turnaround  for  new  technology  [1]. 

One  of  the  greatest  hurdles  to  replacing  all  hydraulic  actuators  with  EMAs  is  the  removal  of  heat 
generated  within  them,  a  consequence  of  the  absence  of  cooling  hydraulic  fluid.  Accurately 
quantifying  the  heat  generated  is  complicated  by  the  highly  transient  and  localized  nature  of  the 
power  demands  of  an  EMA’s  motor,  an  especially  significant  issue  in  aircraft.  Thus,  accurate 
modeling  must  be  dynamic. 
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The  parameters  that  characterize  a  motor,  such  as  resistance  and  inductance,  are  actually 
nonlinear  functions.  In  steady-state  analysis,  these  nonlinearities  can  simply  be  averaged  out, 
and  often  the  parameters  are  considered  to  be  constant  [2],  But  during  transient  analysis  these 
nonlinearities  can  become  critical  [3].  Consequently,  time-averaged,  linear  models  [4]  are 
inadequate  for  thennal  analysis  and  management  designs.  So,  proper  treatment  of  the  problem 
requires  dynamic,  nonlinear  analysis. 

The  motor  parameters  can  be  detennined  by  careful  experimentation  [3-4],  but  this  is  a  time 
intensive  process  and  can  be  expensive  especially  when  considering  multiple  potential  designs 
for  a  motor.  Jens  Otto  suggested  using  a  reduced  order  coenergy  model  [5],  which  yielded 
results  comparable  to  FEM  results.  Likewise,  Roshen  considered  more  involved  empirical 
formulas  which  accounted  for  excess  eddy  current  losses  [6]. 

Through  FEM  modeling,  the  nonlinear  parameters,  such  as  self  and  mutual  inductances  of  the 
three  phase  windings,  can  be  detennined  [7-8],  Though  dynamic  FEM  can  be  very  accurate  [5], 
such  detailed  simulation  is  very  slow.  Fortunately,  FEM  is  only  necessary  to  quantify  the 
parameters;  it  is  not  necessary  for  simulating  the  motor  over  lengthy  mission  profiles.  Instead,  a 
nonlinear,  lumped  element  model  (NL-LEM)  can  use  the  parameters  from  the  FEM  model  and 
then  dynamically  simulate  the  motor  just  as  accurately  as  the  FEM  model  but  at  a  much  lower 
computational  cost.  Our  paper  addresses  the  nonlinear  dynamic  modeling  of  a  pennanent 
magnet  motor  and  describes  both  the  control  and  thermal  perfonnance  of  the  motor  in  following 
highly  transient  input  mission  profiles. 

MODEL 

Two  inputs  are  given:  a  time-history  profile  (a  sequence  of  data  points  versus  time)  of  desired 
positions  (stroke  reference  profile)  and  another  time -history  profile  of  load  forces  (force  profile). 
Both  of  these  profiles  can  be  translated  to  the  perspective  of  the  motor  through  the  conversion 
factor  of  the  drive  train  ratio  so  that  we  get  a  rotor  angle  profile  and  a  load  torque  profile.  Then, 
for  each  time  step,  the  appropriate  applied  voltages  to  make  the  motor  follow  the  stroke  profile  in 
the  presence  of  the  load  torque  profile  are  calculated  by  the  control  algorithm.  These  voltages 
are  then  applied  to  the  motor,  at  which  point  the  simulation  enters  a  convergence  loop  to 
determine  what  the  correct  currents  and  motor-generated  torques  would  be  in  the  presence  of 
such  applied  voltages.  With  the  motor-generated  torque  determined,  the  new  rotor  position  can 
be  calculated  for  that  time  step.  Thus,  the  motion  of  the  motor  is  simulated.  Now,  for  each  time 
step,  once  the  currents  are  calculated,  the  power  loss  in  the  windings  can  also  be  found.  These 
copper  power  losses  then  become  the  heat  sources  for  the  thermal  algorithm,  where  temperatures 
for  various  points  (nodes)  in  the  motor  are  calculated.  The  outputs  from  the  simulation  are  the 
actual  stroke  profile  of  the  motor  (as  opposed  to  what  was  desired:  the  reference  profile),  the 
electrical  power  flowing  in  and  out  (regen)  of  the  motor,  and  the  temperature  profiles  for  the 
various  nodes. 

SIMULATION  EQUATIONS 

The  motor  dynamics  are  modeled  by  four  primary  dynamical  equations  in  direct-quadrature 
reference  frame  (dqO)  [9].  The  first  two  of  these  four  are  shown  here: 
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d/ l. 

dt  -OL* 

(1) 

dA 

+ 

&  * 

+ 

sy 

(2) 

where  ud  is  input  direct  voltage,  u  is  input  quadrature  voltage,  id  is  direct  current,  i  is 
quadrature  current,  Ad  is  direct  flux  linkage,  A,  is  quadrature  flux  linkage,  Rs  is  phase 
resistance,  come  is  the  mechanical  frequency  multiplied  by  the  number  of  pole  pairs  pi  2,  and  p 
is  the  number  of  poles.  The  two  variables  Ad  and  A,  can  be  expressed  by 

\  ~  Ld  (f ,  id  X/  +  Apm  (3 ) 


where  XPM  is  the  flux  linkage  from  the  permanent  magnet,  which  is  a  constant.  Therefore,  the 
time-derivatives  of  these  two  flux  linkages  would  be  the  following: 


dAid  _  /  .  \did  dLd{iq,id)diq  .  dLd[iq,id)did  , 
dt  d  'r  d  dt  diq  dt  d  did  dt  d 

d\  _  /.  .  \diq  dLq(iq,id)diq  ,  dLq(iq,id)did  . 

dt  q  q  ’  dt  diq  dt  q  did  dt  q  ’ 


(5) 

(6) 


For  PM  motors  in  EMA  applications,  we  can  design  the  controller  to  keep  the  direct  current 
close  to  zero.  Also,  the  change  in  Ld  and  Lq  will  change  only  slightly  with  current,  for  our 

particular  motor.  Therefore,  in  equations  (5)  and  (6),  the  second  and  third  terms  can  be  dropped, 
leaving  us  with 


dAd 

dt 

d\ 

dt 


(7) 

(8) 


With  these  changes,  equations  (1)  and  (2)  can  be  modified  to  the  form  of  (9)  and  (10),  which 
along  with  equations  (11)  and  (12)  are  our  four  primary  dynamical  equations. 


dh 

dt 


7“k/  ~RJd  +(>JmiLqiq) 


(9) 
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where  rL  is  the  load  torque,  9me  is  the  mechanical  angle  multiplied  by  the  number  of  pole  pairs 
p  12,  I  is  the  rotor’s  moment  of  inertia,  c  is  the  rotor's  coefficient  of  friction  from  windage  and 
bearings,  and  zM  is  the  motor  torque  generated  by  the  magnetic  fields  and  is  given  by 

=  iq\^P  •  Wm  +  h  ■  (Ld  -  A,)]j  •  O3) 

The  motor  losses  are  related  to  the  motor  parameters,  such  as  Rs ,  c  ,  Ld,  and  L  .  Since  losses 

have  an  effect  on  motor  behavior,  when  possible  they  should  be  modeled  through  the  motor 
parameters,  not  calculated  in  post-processing,  as  is  often  done  [6].  The  copper  loss  can  be 
calculated  as 

Pc=(c+h+c  k.  <14) 

or  equivalently 

p„=  <15) 

for  balanced  3-phases.  The  electrical  phase  resistance  is  itself  a  function  of  temperature.  The 
relationship  between  Rs  and  Rs0 ,  the  phase  resistance  at  room  temperature,  is  as  follows: 

(16) 

A 

where  p0  is  the  resistivity  of  copper  at  room  temperature  and  p  is  the  resistivity  of  copper  at 
any  temperature,  T,  given  as 

p  [nQ  •  m]  =  15.4[l  +  0.0045  \{T  -  273  [k])].  (17) 

In  either  case,  note  that  the  resistance  is  the  parameter  directly  associated  with  the  power  loss  in 
the  windings  (the  copper  loss).  Likewise,  windage  and  bearing  losses  are  directly  associated 
with  the  coefficient  of  friction,  c  .  All  iron  losses  (hysteresis,  classical  eddy,  and  excess  eddy 
losses)  could  be  associated  with  nonlinear,  hysteretic  inductances,  Ld  and  Lq . 
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Much  of  the  motor  modeling  research  done  in  the  past  [2]  and  [10]  and  even  recent  work  [11] 
uses  constant  parameter  values  in  simulations,  although  some  did  use  nonlinear  inductances  [3] 
and  [8],  and  [12].  Since  the  BH  curve  is  nonlinear  and  H  is  a  function  of  current,  inductance  is 
also  a  nonlinear  function  of  current.  For  our  model,  we  used  FEM  to  obtain  Ld  and  L  as 

functions  of  /  .  At  this  point  we  are  neglecting  changes  in  id  since  our  control  algorithm  already 

seeks  to  maintain  a  zero  value  for  id .  The  inductance  curves  were  calculated  in  ANSYS  using 

the  magnetization  curve  for  the  steel-1018  soft-magnetic  material  and  the  geometry  of  our  motor. 
The  analysis  was  done  for  zero  direct  current  and  for  a  wide  range  of  quadrature  currents  (Fig. 

1). 


Inductance  as  a  Function  of  Current 


Fig.  1 .  Direct  and  quadrature  inductances  as  functions  of  quadrature  current 
TIME  STEPPING 

The  size  of  variable  arrays  of  data  points  affects  the  speed  of  a  simulation.  Using  prerecorded 
input  profiles  (time-history  arrays  of  stroke  and  load  torque),  very  small  time  stepping  is 
necessary  to  accurately  capture  the  transient  behavior  in  the  motor.  A  uniformly  applied,  very 
small  time  stepping  would  result  in  very  large  data  arrays.  To  keep  our  array  sizes  small  while 
maintaining  small  time  stepping  where  needed,  we  used  a  line  simplification  algorithm  to  reduce 
the  number  of  points  where  little  transient  activity  would  occur.  This  data  compression 
technique  allowed  us  to  reduce  the  input  profiles  to  less  than  4%  of  their  original  data  size. 

That  said,  the  electrical  time  constant  is  often  much  smaller  than  the  smallest  time  step  of  the 
input  profiles.  Since  the  frequency  of  the  electrical  component  of  the  simulation  is  directly 
related  to  the  first  derivative  of  the  input  angular  position  profile  for  the  motor,  we  can 
dynamically  scale  the  electrical  time  step  non-unifonnly  [13].  Therefore,  we  had  two  loops,  one 
macro  loop  that  stepped  through  the  elements  of  the  input  stroke  profile  and  a  nested  micro  loop 
with  very  small,  but  variable  time  steps  for  the  electrical  equations.  The  result  of  this  apparent 
complication  was  a  significant  speed  improvement. 
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FIELD  ORIENTED  CONTROL  (FOC) 


A  thorough  analysis  of  motor  heat  generation  should  include  how  the  motor  is  driven.  Our  field 
oriented  control  is  shown  in  Fig.  2,  where  the  numbers  in  parenthesis  refer  to  the  equations 
shown  further  below.  The  control  is  an  adaptive  closed-loop  control,  which  is  a  type  of  feedback 
linearization  of  the  speed  and  position  errors. 


Fig.  2.  Field  oriented  control  (FOC)  of  PM  motor 

Our  control  equations  are  derived  from  the  primary  simulation  equations  (9)  through  (12). 
Starting  with  equation  (1 1)  and  noting  that  in  (13)  i  can  be  factored  out,  we  get 

\  —  =  ^L+-(I^me+C(°lne)’  (  1  8) 

C  P 

where,  the  angular  acceleration  multiplied  by  the  number  of  pole  pairs  p/2,  ame ,  is  defined  as 
time  derivative  of  a>me .  Equation  ( 1 8)  can  be  rearranged  to  be 


C  = 


TL+-{lame+ccome) 
P 


(19) 


At  this  point  we  reinterpret  ame  which  is  dcome  /  dt  to  be  ef0  /At ,  where  At  is  a  finite  time  step 
and  e(0  is  the  error  in  0)me,  with  a  tuning  coefficient,  kp  ,  in  front.  This  gives 
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tl  + 


Ikp  r. +  cco »* 

At 


z„  = 


'  M 


(20) 


9]e  s$c 

where  iq  is  the  desired  current  for  the  next  time  step  and  come  is  the  desired  next-step  angular 
speed  equal  to  came  +  kpe0)  .  Using  (1 1)  to  substitute  for  r,  in  (20),  we  arrive  at  (21): 


h  = 


2 

P 


rM - Iam„ - cco„ 


2 

P 


+  ■ 


Ik„ 


At 


+  cco„ 


(21) 


Recalling  that  [come*  -  come )  equals  kpe0) . ,  (21)  simplifies  to 


* 


(22) 


This  equation  will  track  the  desired  motor  velocity  well,  but  there  will  be  some  rotor  angle  drift 
over  time.  To  counter  this  drift  error,  we  added  a  theta  error  term,  Uty  /  At : 


* 


l9  + 


'  At 


+ 


I 

(  .  ea,me 

k  „ - 

-  -  a . 

+  c\kne„  ) 

^  p  At 

me  j 

V  p  come  J 

.(23) 


The  kp  and  k  coefficients  are  the  only  PI  (proportional  integral)  values  that  need  to  be  tuned  in 
the  control  algorithm.  The  coefficient  kp  scales  the  controller's  reaction  to  the  error  in  0)me ,  and 

kt  scales  the  controller’s  reaction  to  error  in  0me ,  the  motor  angle  multiplied  by  the  number  of 

pole  pairs.  Together,  these  two  coefficients  help  us  set  how  quickly  the  controller  will  try  to 
eliminate  the  particular  error. 
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For  maximum  efficiency,  it  is  generally  desired  that  the  direct  current  be  zero  since  any  flux 
generated  in  the  direct  axis  would  not  contribute  to  torque.  However,  direct  current  does  not 
necessarily  directly  translate  to  direct  flux.  Direct  current  can  contribute  to  useful  torque  if  the 
direct  and  quadrature  inductances  are  not  equal,  as  given  in  (13). 

$ 

Most  of  the  time,  however,  the  desired  direct  current  set  point,  id  ,  should  be  zero.  With  the 
desired  currents  detennined,  the  desired  voltages  can  be  found: 

*  *  *  di .  *  *  *  _  . 

C/  =  Rsld  +  Ld  ~Jt - <*>me  Lq  \  (24) 


U„  =RJ,  +L„ 


di„ 


dt 


*  *  .  *  * 

'  mme  Rd  id  ® me  PM  > 


(25) 


where  Ld  and  Lq  are  based  on  the  desired  current  values.  It  should  be  noted  that  while  many 
of  the  variables  are  updated  based  upon  the  desired  currents,  rM  was  not  updated  for  the  control 
equations.  This  is  a  reasonable  simplification  since  for  a  round  rotor  XPM  is  independent  of  6me , 
id  should  normally  be  zero,  and  the  i  factor  is  divided  out: 


—  =  (26) 
h  4 

For  the  non-round  rotor  case,  a  direct  interpolation  method  with  a  lookup  table  could  be  used  to 
include  the  dependence  on  theta. 

Even  after  the  desired  control  values  were  calculated,  the  reality  of  limits  (current  and  voltage) 
needed  to  be  included.  Instead  of  directly  capping  those  values,  we  realized  that  the  system 
converged  more  quickly  and  ran  more  smoothly  if  the  functions  used  were  all  smooth  [14].  For 
our  case,  we  used  a  smooth,  segmented  curve  of  first  and  second-order  polynomials  to  define  the 
saturation  of  voltages. 

TEMPERATURE  CALCULATION 

Using  a  lumped-element  thermal  model,  the  temperatures  at  12  different  locations  (Refer  to 
nodes  1,  3-10,  and  12-14  in  Fig.  3)  in  the  motor  were  calculated.  Similar  work  has  been  done 
before  [15],  though  with  only  one  node.  The  source  of  heat  for  the  temperature  calculations  was 
the  copper  power  loss  in  the  windings.  The  temperature  distribution  of  the  motor  was  modeled 
using  a  lumped-element  network  (Fig.  3),  which  is  limited  to  what  is  inside  the  motor  assembly 
(indicated  by  the  dashed  line  in  Fig.  3).  The  heat  transfer  coefficient  used,  h  ,  was  1.45  W/K  m2. 
The  values  of  the  thermal  resistances  and  capacitances  were  tuned  by  comparing  the  results  of 
the  lumped-element  model  with  those  of  an  ANSYS  FEM  model  of  the  motor. 


Detailed  Lumped  node  model  Of  EMA  (transient) 


I.  copper 

3.  stator 

4.  case  and  back  cover 

5.  case  surface 

6.  shaft  body 

7.  shaft  front  end 

8.  shaft  back  end 

9.  magnet 

10.  air  gap 

II.  connection  gear 

12.  front  bearing 

13.  back  bearing 

14.  front  cover 

15.  moving  mechanism 

16.  wing  surface 

17.  gear  box  case 


Motor  assembly 


Fig.  3.  Thermal  lumped  element  model 


This  model  is  generic  with  multiple  heat  sources.  However,  for  our  simulations  we  considered 
only  the  copper  loss  which  is  the  predominant  heat  source  in  our  motor,  and  the  other  heat 
sources  were  set  to  zero.  A  few  of  the  locations  of  the  temperature  nodes  are  shown  in  Fig.  4. 
Node  5  is  the  case  surface  and  connects  directly  to  the  ambient  air  which  serves  as  the  heat  sink. 


Airgap  (Node  10) 
Magnet  (Node  9) 


Rotor  core 
(Node  6) 


Case 
(Node  4) 

Case  surface 
(Node  5) 


Stator  iron 
(Node  3) 

Copper  winding  (Node  1) 


Fig.  4.  Thermal  node  locations  on  motor 

RESULTS 


The  input  (desired)  stroke  profile  and  the  motor's  actual  stroke  profile  are  shown  in  Fig.  5.  These 
profiles  are  made  of  18  repetitions  of  a  5  minute  profile. 
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Mechanical  Stroke  Following 


"0  10  20  30  40  50  60  70  80  90 

Time,  f  [min] 

Fig.  5.  Desired  (blue)  and  actual  (green)  stroke  profdes 

The  motor  followed  the  desired  stroke  profile  very  well  which  is  why  barely  any  blue  can  be 
seen.  The  input  load  torque  profile  is  shown  in  Fig.  6. 


Load  Torque 


‘20  10  20  30  40  50  60  70  80  90 


2, 

1.5 

1 


u  0 
0) 


Time,  t  [min] 

Fig.  6.  Input  load  torque  profile 

We  have  found  that  the  control  design  influences  the  amount  of  heat  generated,  especially  the 
peak  power  lost  to  heat.  When  k  was  set  to  0.05  s  ,  the  mean  power  loss  in  the  windings  for 

the  given  input  profile  was  about  0.99  W,  the  peak  power  loss  was  12.05  W,  and  the  error  (a 
normalized  sum  of  differences)  in  stroke  (motor  position)  following  was  0.0817%.  In  Fig.  7,  we 
zoomed  in  on  the  stroke  following  to  show  6  seconds.  Here,  a  little  overshoot  can  be  seen.  We 
will  call  this  performance  loose. 
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Mechanical  Stroke  Following 


Time,  t  [min] 


Fig.  7.  Zoomed-in,  6-second  view  of  loose  stroke  following  (k  =  0.05),  desired  (blue)  and 
actual  (green) 

When  the  k  value  was  increased  to  0.7  s'1,  the  mean  power  loss  increased  to  about  1  W,  the 

peak  power  loss  rose  to  20.65  W,  but  the  stroke-following  error  dropped  to  0.0248%  (Fig.  8), 
meaning  tighter  control. 


Mechanical  Stroke  Following 


Time,  t  [min] 


Fig.  8.  Zoomed-in,  6-second  view  of  tight  stroke  following  ( kp  =  0.7),  desired  (blue)  and  actual 
(green) 

Although  the  mean  power  loss  for  these  two  simulations  was  about  the  same,  there  was  a  large 
difference  in  the  peak  power  loss.  In  our  simulations,  we  found  that  a  kp  value  set  too  high  can 

cause  additional  power  loss.  The  additional  power  loss  comes  from  the  motor  working  too  hard 
to  follow  the  input  profile  exactly. 

The  copper  power  loss  (for  kp  =0.1  s'1)  in  the  motor  for  a  90  minute  run  is  shown  in  Fig.  9. 
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Power  Loss  in  Windings 


Time,  t  [min] 

Fig.  9.  Power  loss  in  windings 

The  mean  power  loss  was  0.98  W  and  the  peak  power  loss  was  1 1 .7 1  W.  It  is  worth  noting  that 
the  peak  power  loss  was  12  times  higher  than  the  mean  power  loss.  This  points  to  the  necessity 
of  an  accurate,  transient  analysis  versus  a  steady-state  analysis  since  such  surges  of  power  can  be 
the  cause  of  motor  failure. 

Of  great  interest  to  power  systems  engineers  are  the  electrical  input  power  and  the  regen  power. 
Our  simulation  outputs  both  of  these  quantities  (Fig.  10). 


Electrical  Power 


Fig.  10.  Input  electrical  power  (blue,  positive)  and  regen  power  (red,  negative) 

The  temperature  profiles  for  various  nodes  are  shown  in  Fig.  1 1 .  The  highest  temperatures 
pertain  to  the  nodes  closest  to  the  copper  windings.  The  ambient  temperature  was  held  at  22  °C. 
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Fig.  11.  Temperature  profiles 

Though  this  motor  was  put  under  a  fairly  light  load,  the  critical  aspect  of  nonlinear,  transient 
simulation  is  being  demonstrated.  The  3  °C  temperature  rise  can  be  simply  explained  by  energy 
balance: 


Q  =  ATmcp  (27) 


where  the  total  energy  lost  to  heat,  Q ,  was 
Q  =  0.98  W  (60s-  90)  =  5292  J  (28) 

The  total  temperature  rise,  AT ,  was  3  °C;  the  mass  of  the  motor,  m  ,  was  4.4  kg;  and  the  specific 
heat  capacity,  cp ,  was  estimated  as  400  J/kg  °C.  Therefore, 


ATmcp 


(3°CX4.4kgj  400 


x 


kg°C 


= 5280 J 


(29) 


which  is  only  a  0.23%  difference  from  (28). 

Although  the  speed  of  the  computer  will  certainly  affect  the  run-time  ratio  (real  time/simulation 
time)  as  will  the  particular  profile  being  simulated  since  highly  transient  portions  of  the  input 
profiles  require  more  computations,  on  the  Gateway  Centrino  Duo  2.16  GHz  laptop  we  used  and 
with  our  inputs,  the  real-time  ratio  was  over  20. 

CONCLUSIONS 


FEM  is  a  powerful  tool  for  characterizing  motor  parameters,  and  the  combination  of  FEM  and 
NL-LEMs  provides  for  accurate,  fast  simulations  of  EMAs  under  transient  conditions.  This 
paper  shows  our  integration  of  dynamic,  nonlinear,  lumped-element  modeling,  tuned  by  FEM, 
with  field  oriented  control,  copper  loss  and  regenerative  power  calculation,  and  multi-node 
thermal  modeling  of  the  transfer  of  generated  heat.  Since  our  model  runs  so  fast,  it  could  be  used 
in  hardware-in-the-loop  simulation.  We  detailed  the  primary  dynamical  equations  of  the  motor, 
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and  showed  how  these  are  used  to  simulate  the  motor  and  to  build  the  control  for  the  motor.  We 
showed  how  the  thermal,  lumped-element  network  was  laid  out  and  where  several  of  the  network 
nodes  were  placed  in  the  motor.  And  we  showed  how  control  can  have  an  effect  on  the  peak 
power  loss  which  dissipates  as  heat  and  on  the  stroke-following  performance.  Finally,  we 
showed  how  the  inputs  of  desired  stroke  profde  and  load  torque  profde  led  to  the  outputs  of 
actual  motor  stroke  profile,  input  power  and  regenerative  power  profiles,  and  temperature 
profiles. 
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NOTATION 
VARIABLES 
u  voltage  [V] 

1  current  [A] 

R  electrical  resistance  [Q] 

L  inductance  [H] 

2  flux  linkage  [Wb] 

6  angle  [rad] 

co  angular  frequency  [rad/s] 

a  angular  acceleration  [rad/s  ] 

r  torque  [N  m] 

p  number  of  poles 

/  rotor's  moment  of  inertia  [N  m  s“/rad] 

c  rotor's  coefficient  of  friction  from  windage  and  bearings  [N  m  s/rad] 

P  power  [W] 

e  error,  [rad/s]  or  [rad/s  ]  depending  on  the  units  of  the  subscript 

k  scaling  coefficient  for  error,  [1/s]  for  kp  and  [A  s/rad]  for  kt 
h  heat  transfer  coefficient  [W/K  m  ] 
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dt  infinitesimal  time  step  [s] 

At  finite  time  step  [s] 

SUBSCRIPTS 

d  direct  axis 

q  quadrature  axis 

s  stator 

me  mechanical  scale  multiplied  by  number  of  pole  pairs 
PM  permanent  magnet 
M  motor 

L  load 

a  phase  a 

b  phase  b 

c  phase  c 

cu  copper 
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